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Abstract 

Using the notion of aggregation work, we construct a system of differential 
equations for the aggregation number of micelles which is a function of the pa- 
rameters of micellization {parametric equations). There are explicit solutions 
for two important models of spherical micelles. Based on these solutions, we 
obtain an analytical expression for the equilibrium concentration of surfactant 
monomers and consequently for the whole spectrum of equilibrium concentra- 
tions of molecular aggregates in this framework. Accuracy of these expressions 
is discussed, and they are applied on an example of micelles formed by sodium 
dodecyl sulfate. 



Introduction 

The process of micellization is an important and interesting problem, whose 
mechanism is not yet understood. Its complexity does not allow to describe it 
completely. If the concentration ni of surfactant monomers (amphiphiles) ex- 
ceeds the critical micellization concentration, the monomers can form aggregates 
called micelles. There are different structures of micelles (spherical, cylindrical, 
disk-like, inverse, etc., see [|l|), and they are classified by ni. For a relatively 
small concentration ni, the spherical structure is the most favorabl^ There 
are two essential advantages in this case. First, the concentration of monomers 
is small and therefore one can apply some useful approximations. Second, the 
geometry of micelles is described by the only parameter - the radius of the 
sphere. Moreover, this radius can be expressed in terms of aggregation num- 
ber (the number of monomers which form an aggregate). Consequently, it is 
natural to use methods of the nucleation theory and, particularly, the notion of 
aggregation work. One can find its description and applications in that 
investigate thermodynamic properties of micelles. 

We propose a method which allows to determine the dependence of the 
equilibrium concentration of monomers and micelles on the essential parame- 
ters defining the process of micellization. We are going to establish parametric 

^ E-mail: Donis.Grebenkov@polytechnique.fr 

^ Here we will not consider the supplemental conditions which can be imposed by packing 
constraints, see 
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equations which can be solved exphcitly for the different model expressions of 
the aggregation work. Therefore, we will able to determine equilibrium concen- 
tration of monomers as a function of external parameters (e.g., temperature) in 
the framework of a given model. Consequently, it will allow to relate different 
models of spherical micelles among themselves and to compare them with regard 
to the experimental data. 

The first part is devoted to introduction of the basic notions and the essential 
assumptions which will be used to establish parametric equations in the second 
part. The third and the fourth parts give solutions of these equations for the 
drop model and Grinin's model of spherical micelles. The final part describes 
some consequences and the possible generalizations of this formalism. 



1 Basic notions and essential assumptions 

Nucleation theory is based upon the notion of aggregation work F^,, i.e. minimal 
work required to form an aggregate from v surfactant monomers (amphiphilcs), 
where v is referred to aggregation number. Function F^, essentially defines the 
dynamics of the system, and completely determines equilibrium distribution of 
aggregates by virtue of the aggregation number according to Boltzmann's 
law, 

= nie'^"", (1) 

where ni is the concentration of monomers]^ Note that Fi = because the 
monomers already exist. Function F^, depends on the aggregation number v, 
concentration of monomers in solution ni , and parameters defining the micelliza- 
tion (e.g. temperature, molecular properties of solvent, structure of monomers, 
etc) which are denoted as ai...aK- Note that the aggregation work F^ does not 
depend on the concentrations (j^ > !)■ Indeed, we are working with spherical 
micelles, thus the concentration of monomers is relatively small and highly dom- 
inates the concentrations of dimers, trimers, etc. Hence, we can assume that 
aggregates are formed due to gradual attachment/detachment of monomers. In 
other words, we neglect the possibilities of attachment/detachment of dimers, 
trimers, etc. to the aggregate. Moreover, for the small concentration (i.e. the 
dilute solution) the dependence of function F^ on ni is simplyQ 

^ We shall use dimensionless quantities. Particularly, all the concentrations are expressed 
in convenient units: for example, in units of the critical micellization concentration; the 
aggregation work is measured in units of kT, where k is Boltzmann constant, T is absolute 
temperature. 

This is a direct consequence of the differential form of as difference of chemical 
potentials 

dF„ = {^y — iJ.o)di', 

where /i^ is the chemical potential of any surfactant "inside" the aggregate consisting of v 
molecules. Obviously, is independent of the concentration of monomers in the solution 
and depends only on the internal structure of aggregate, /^o is the chemical potential of any 
surfactant in the solution. For the dilute solutions there is a simple expression for /^o 

/iO =lnni+x(T), 

where xC^) ^ certain function of temperature (see |7|). 
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=G, - (j.- l)lnni, (2) 

where new function Gi, does not depend at all on the concentrations and is 
defined only by the parameters ai..aK and v. The condition Fi = implies 
Gi =0. 

We make an important assumption of the closure of the system, i.e. the 
total number of amphiphiles (in monomers and aggregates) N is fixed. This 
is an essential condition because one can imagine some processes where this 
number changes. From now on, N will be a fixed external parameter. In this 
case, the law of mass conservation can be written as 

oo 

N^^^i^n,. (3) 

Note that in real physical problems the summation is usually limited by a certain 
value I'm- However, it is convenient to use this notation, remembering that the 
finite sum can be "complemented by zeros" (which is equivalent to rapid growth 
of Fi, with increase of ly) . For the convergence of series it is sufficient that F^, 
tends to infinity when ^ cx), or that Gi, increases faster than a linear function 

of I/. 

As the mechanism of micelle formation is very complex, it is reasonable to 
use model expressions for G^, which should represent the essential properties of 
micellization. We represent function as a sum of the fractional power^ of 

M 

G. = ^ c^i.'"", (4) 

m=0 

where p > is a constant for the particular modelQ, the coefficients Cm of this 
decomposition depend on the parameters ai...aK, and there are only positive 
powers in the sum^. Note that function Gi, has a physical sense only for natural 
numbers ly. However, the formulae (^ allows us to interpolate G^ for any real 
number v > 1. Obviously, one can imagine different interpolations but we 
shall use the simple one|^ given by (Q). Furthermore, we shall differentiate some 
quantities by etc. 

The important statement of the nucleation theory is following: every mini- 
mum of aggregation work corresponds to a stable (or quasi-stable) state of the 

^ This representation is quite reasonable: for example, in the problem of condensation of 
vapor the aggregation work is the sum of the powers of u^^'^' , in Grinin's model of spherical 
micelles — the sum of the powers of u^^^. 

^ As we use the model expressions, we can think of p as a rational number. If the sum 
includes the linear term, we can take p = 1/mo with a natural number mo- 

^ The negative powers of u in the decomposition have no practical use because we are 
interested in the quantity for large u. 

* Nevertheless, sometimes it is convenient to change a little this canonical interpolation. 

Indeed, the limit of Gi, when v — > 1 is equal to c = X/m— 0^™' ^'^'^ nothing demands that 

it equals to zero (for example, see the drop model of spherical micelles, part 3). Therefore, 

we can change the expression Mj on the interval (1, 2) to obtain lim Gi, = 0. The simple 

1/ — ^1 

way is to add function —c/u^ to Gv where x is a large positive number. We are interested in 
function Gi, for large u thus the addition —c/u^ has no any influence on our results. 
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system of aggregates. We are interested in the formation of micelles, therefore 
there should be a "deep" minimum of F^, (potential well). Moreover, this min- 
imum is assumed to be unique and it will be denoted as i^s- Indeed, we are 
working with a small concentration of monomers, and consequently, there is 
only one stable state - spherical micelles. Using the relation (||), one can write 
the equation on Vs as 

= ^^(i^s) = ^^(I's) - Inni. (5) 

Physically, the quantity Us is the aggregation number of micelles, Vg ^ \. In the 
practically important case the number of amphiphiles in micelles is comparable 
(or more) to the concentration of monomers. 

Using all previous assumptions we can establish our formalism which defines 
the dependence of equilibrium concentration of monomers on the parameters of 
micellizaton. 



2 Parametric equations of the theory of forma- 
tion of spherical micelles 

Let us introduce a set of functions fa{x) assuming x > 0, 

oo oo 

/„(a;) = ^j."e-«"a;'' = ^e-''(''), (6) 



where new function hiy) is 



hiy) — Gi, — vlnx ~ alnh'. (7) 
There are three important properties : 

1) All functions fa{x) are defined for any positive x because of a rapid growth 
of function h{i') (since increases faster than v). Moreover, fa{x) > is true 
for any a: > 0. 

2) Simple verification shows that 

X—g^=fa+l[X). (8) 

Hereof it follows that all fa{x) are smooth monotonously growing functions, 
with derivatives of all orders. 

3) For a i^s function h{v) has a minimum which is with high accuracy]^ 
equal to Vg. 

Let us represent fa{,x) as a sum of three terms which correspond to mono- 
mers, small aggregates and micelles, 

1^0 oo 

(x) = X + ^ -1- ^ e-^^") = X- + /o (x) + /„ (x) . (9) 



^ The equation for the minimum of h{u) is derived by subtracting a/u in the right-hand 
side of the equation (Bh. Under condition a ^ Us this term can be neglected by comparison 
with the positive powers of u. 
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It is necessary to say that the choice of this decomposition is artificial: i^q 
depends on the given model. Usually vq can be chosen as the maximum of F^, 
i.e. a point where rii, has its minimum. As we take interest in formation of 
micelles, it is reasonable to suppose that the contribution of micelles dominates 
over the contributions of small aggregates, 

(10) 

We neglect functions /^(x), and the sum (^) becomes 

Ux)^x + f^{x). (11) 

Now we calculate an approximation to function fa{x) (which corresponds to 
micelles) for small values of a. For this purpose we represent the corresponding 
sum in (|^) as an integral, decompose hiv) in the vicinity of its minimum Vg up 
to the square terms, and calculate the obtained gaussian integral after replacing 
the lower integral limit by —oo, 

/„(a:)«exp[-/3(x)], P{x)^h{v,) + Unh"{vs)~\H2^). (12) 

Due to the independence of Vs of a, we obtain an important relation between 
functions /„, 



whence we conclude 



h = {^sT'^h- (13) 

In terms of functions /« (x) the law of mass conservation (|^) can be written as 
equation for rti, 

h{x)^N. (14) 

Due to the continuity and monotone growth of fa {x) , for any A'^ there is unique 
solution of (0) which is equal to ni, 

h{ni)=N. (15) 

Substituting (|lj) and ( [l5| ) in (|Tl|), we obtain the transcendental equation for 
the equilibrium concentration of monomers ni (cf. [^), 

N = ni +exp[-/3(ni)]. 

Practically this equation leads to the same difficulties in solution as the initial 
equation (p^). Therefore we continue to develop our method. 

Let us consider the equation (^4|). Once this equation is solved, we shall 
find the dependence of x (i.e. ni) on the parameters ai...a/f. One can treat 
this dependence as a trajectory in the phase space based on parameters ai...ax, 
where function fi{x) is "integral of motion" (i.e. it is a constant on this trajec- 
tory). Fixing all parameters except one and differentiating (|l^ ) by this single 
parameter a^, we express the derivatives dx/duk according to 
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dx 
dak 



dfi 



dx 



-1 



oak 



The derivative in denominator can be found with the help of (^. Now we 
calculate the derivative in nominator, 

9fi ^ -G. .f dG^ 
ve "x 



dak V ^o-k 

Taking into account the representation of Gv as a sum of fractional powers of 
one has 

^ = Vc™,,.^"™, c™.. = ^. (17) 
oak ^ oak 

in 

Collecting the previous expressions together, dividing by x and bringing it under 
the sign of differential, we obtain 

= J / \ / . Cm,kjpm+1- (18) 

dak hyx) ^ 

Now we represent functions fa according to (|^), where all the terms fa are 
expressed by 



1^ ) and all the terms /°(x) are neglected by virtue of (|i(]|), 

d\nx 1 
dak x + vsfi 



x^^c^.k + 



ra / 



In the first part we demanded that the number of amphiphiles in micelles (/i) 
has the same or much greater order than in monomers (r^i), therefore we can 
neglect]^ x and x^^Cm,k by comparison with Vsfi and /i Z^m ^m.fcJ^f" in 
nominator and denominator respectively. We have 

m 

With equation (^ we obtain a closed system of equations for two unknown 
functions : In x and Vg ■ Differentiating the equation (||) by ak , 



din; 
dak 



we establish one differential equation for lyg. After simplification, this equation 
gets a form 



, I]Cm,fc(/5m - 1) 



dak J2 Cmpm{pm - l)z^r" ^ ' 



(fc = I..K). (20) 



So, we have obtained a set of equations which define the dependence of the 
aggregation number of micelles Vg on the parameters of the theory. We shall 

^'^ Here we make second approximation, after the gaussian one in calculation of fa- The 
essential simplification brought by this operation is the independence of derivatives d\nx/dak 
of /l. 
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call these equations parametric. Its solution is connected with the equilibrium 
concentration of monomers ni via (^). 

If we are interested only in dependence on one specific parameter Ofc, it is 
sufficient to solve one of the equations ( po| ) for this parameter. However, in 
this case there appears an arbitrary function of the other parameters. On the 
contrary, if we search for a complete solution, it is necessary to solve all the 



equations (20) and to accommodate all appearing arbitrary functions. Note 
that in the case of a large number of parameters this algorithm can be very 
complicated. 

We make another important remark. If one adds an arbitrary function 5(1^) 
to the expression of Gi, which does not depend on the set of parameters ai...aK 
and does not break a convergence of the series (j^), then the equation ( poj ) does 
not change. Indeed, function g(v) disappears after differentiation by the pa- 
rameters ai-MK- Consequently, g{v) has no influencef^ on a solution of the 
parametric equations (pO|). On the contrary, the addition of g{v) affects the 
equation (|^) , and the expression of equilibrium concentration of monomers will 
change according to the formula 

n[ ^ niexp[-g'{iys)], (21) 

where ni is a solution of ( pO[ ) without giv). This simple proposition is very useful 
from the practical point of view. We see that the characteristic i>s is universal, 
its equation ( [20[ ) is independent of a small variations of function G^. On the 
contrary, the equilibrium concentrations of monomers and micelles depend on 
the derivative of .g(i^) according to (^l]). 

The common property of the set of equations ( po|) is their independence of 
the linear term v which disappears by virtue of a factor {pm — 1) when m equals 
to mo — l/p corresponding to a linear term (in the case when l/p is not natural 
number, a priori there is no linear term). Nevertheless, sometimes^ there is a 
dependence of the equation ( ^ ) on the coefficient Cmg by virtue of the term cq 
which can contain this coefficient. However, usually one can neglect the term 
Co by comparison with j^f™ {i>s ^ 1). Note that every case should be analyzed 
particularly. 

What is the qualitative information that is contained in the equation (pO|)? 
It seems possible to use this equation (and the whole method) to verify the fol- 
lowing important hypothesis in the theory of micellization. As it was mentioned 
above, the exact expression for the aggregation work is unknown, therefore we 
use the model expressions (see, for example, (^T])). However, it is natural to 
suppose that all physically measured quantities depend on global characteris- 
tics of the aggregation work (such as minimum of , depth of this minimum, 
height of activation barrier, etc), rather than on a choice of its particular ex- 
pression. If this hypothesis works, it is sufficient to take a simple expression for 

satisfying to the qualitative requests and to solve the problems with the use 
of this expression. 



It is more accurate to say that function has no influence on the analytic form of 

solution, but it changes, of course, initial conditions, and consequently, some constants of 
solution. 

See, for example, Grinin's model, part 4. 
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3 Solution of the parametric equations for the 
drop model of spherical micelles. 

Let us consider a drop model of spherical micelles, firstly elaborated by Is- 
raelachvili, Mitchell and Niham (see and reformulated independently by 
Rusanov in terms of the aggregation work (see [^). The basic ideas of this 
model are simple : 

1) hydrophilic electrically charged heads form a surface of micelle (sphere); 

2) hydrophobic tails are very flexible, and they form a liquid drop in the 
micellar core. 

There are three terms in the aggregation work : 

- bulk term (proportional to the number v of amphiphiles in the aggregate) ; 

- surface term corresponding to the attractive hydrophobic forces (and pro- 
portional to sv, where s is the surface area per amphiphile ; sv ~ i^^/"^ for the 
spherical micelles) ; 

- surface term corresponding to the repulsive electrostatic forces (propor- 
tional to v/s v'^l^ as it has been shown by Tanford, JlOt). 

Leaving out the details, we write only the final expression for the aggregation 
work for this model. 



=6ii.4/3„ |^_y^j^^^g^2/3^ (22) 

where the parameters^ h\ and 63 are defined through the physical characteris- 
tics, 

_ {ez)H _ 47rA^7o 

- 8^£oeA2fcT' - ^I^' ^^^^ 

where A = (3i;/47r)^/"^ is radius of a sphere which has the same volume as 
one hydrocarbon chain, ez is electric charge of one amphiphile, Eq is dielectric 
constant, e is dielectric permittivity, 8 is separation of the capacitor planes, 70 
is surface tension of the hydrocarbon-water interface. 



The expression (22) defines the coefficients of decomposition of the aggrega- 



tion work by the powers of v with p = 1/3, 



Co = 0, Ci = 0, C2=fe3, C3 = --v/26163, C4 = 5i. 



In the framework of this model the parametric equations (|2C|) can be solved 
exactly, without any approximation (however, the parametric equations them- 
selves remain approximate ones). Substituting the values of c„j into the equa- 
tions (po|), we obtain 



dv,, 3 vt^^ 



dh, 9 oj, 1/3 L -1/3' 

dv, 3 v^^^ 



1^ We use Rusanov's notations. 



dh^, 9 01, 1/3 v -1/3' 

""3 ^ 26iZys -Osl's 



(24) 
(25) 
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The substitution 



leads to 

^ Lpjip"^ - 1) u ^ ^ 'fii'fi^ ~ ^) 

and we obtain the sohition 



(^2(<^2_1) = 73(&3)6i, (26) 
V^ifi'-l) = li{hi)bf. (27) 

Equaling the right-hand sides, we relate functions 71 (&i) and 73(^3), 

7i(&i)fcr' = 73(?^3)&i 

As the parameters bi and 63 are independent, both sides of this equality are 
constant, whence 

7i(fei)=7Wfci, lM=l{N)b-\ 

where a constant 7(A^) can be determined numerically for any N (see table 1). 
Substituting function 71(61) in the equation (^), one finds the dependence of 
the aggregation number of micelles Vg on the parameters of micellization. 



vs^{ab^/2bif'^, a ^ I + ^ \ + A^{N)bib-\ (28) 

The equilibrium concentration of monomers can be found by substitution of 
28|) in the equation (^), 

^y2MlU/4-a-i/4l'V (29) 



ni = exp 

Now it is useful to evaluate the values of parameters 61 and 63. We can 
obtainf] 

61 w 2, 63 « 30 (30) 

for sodium dodecyl sulfate {C'i2H250S03Na) in water at the ambient temper- 
ature {T = 300i4r) using characteristic values of the following parameters (see 

- volume of the hydrocarbon tail of an amphiphile, v — 0.350 nm^ ; 

- interfacial free energy per unit area of aggregate, 70 — 50 erg/cm^ ; 

- optimal surface area per amphiphile, oq ~ 0.62 nm^ (this value was used 
to calculate 61, see for details [|l|). 

We have determined the dependence of equilibrium concentration of mono- 
mers on the parameters 61 and 63. The dependence on the total concentration 
of amphiphiles N is contained in the constant 7(A'^) which can be numerically 

Note that these estimations are quite enough for our purposes. Indeed, we are working 
with the simple model, therefore we can not expect a high accuracy of predicted results. 
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calculated for any N. Indeed, it is sufficient to find nl for one set of bl, 6° 
and N by solving the equation ^ or , then to substitute the value of Ui in 
the expression (|2^) and to express the constant 7. With the help of numerical 
simulations we obtain the values of 7(A^) for some N represented in the table 1. 



N 


2.0 


2.5 


3.0 


3.5 


4.0 


4.5 


5.0 


ni 


1.9865 


2.1159 


2.1471 


2.1648 


2.1771 


2.1865 


2.1942 


1 


-32.774 


-22.852 


-20.486 


-19.156 


-18.232 


-17.525 


-16.953 



Table 1. The values of constant "f{N) and of equilibrium concentration of 
monomers ni as functions of N . For these simulations we have used bi — 2.0, 
63 — 30 corresponding to sodium dodecyl sulfate in water. 

Note that the aggregation work Fi, for this model grows rapidly for small 
V ; thus dimers, trimers and other small aggregates are completely negligible. 
Consequently, the number of amphiphiles in micelles can be well approximated 
as — rii. 

The essential advantage of the expression ( p9| ) is its explicit form which 
allows to analyze dependencies of all the equilibrium concentrations on param- 
eters &i and &3. Moreover, this formula has a quite high accuracy. The table 2 
represents its relative errors (in per cent) when parameters bi and 63 are widely 
varied. 



&3\6l 


1.0 


1.5 


1.6 


1.7 


1.8 


1.9 


2.0 


2.1 


2.2 


2.3 


2.4 


2.5 


28 


0.9 


0.5 


0.3 


0.1 


0.0 


0.3 


0.5 


0.7 


1.0 


1.3 


1.6 


1.9 


29 


0.9 


0.6 


0.4 


0.3 


0.1 


0.0 


0.2 


0.4 


0.7 


0.9 


1.2 


1.5 


30 


0.9 


0.7 


0.6 


0.4 


0.3 


0.2 





0.2 


0.4 


0.9 


0.8 


1.0 


31 


0.9 


0.7 


0.7 


0.6 


0.5 


0.3 


0.2 


0.1 


0.1 


0.6 


0.4 


0.6 


32 


0.9 


0.8 


0.7 


0.7 


0.6 


0.5 


0.4 


0.3 


0.2 


0.0 


0.1 


0.2 



Table 2. The relative errors (in per cent) of the formula ^{\ ) when bi changes 
from 1.0 to 2.5, 63 changes from 28 to 32. Generally, these errors are less than 
1% (beside of the right upper corner marked in italic that is not physical, see 
below). 

Note that the variation of 63 from 30 to 32 (28) can be caused by decreasing 
(increasing) the temperature by 20°C (i.e. the domain of variation of temper- 
ature is between 5°C and 45°C). On the contrary, the parameter 61 can be 
varied in the wide enough region. Indeed, if one adds salts (for example, NaCl) 
in solution, they will screen electrostatic repulsion of amphiphile heads, i.e. the 
parameter bi will be decreased. Normally, it is difficult to increase 61 consider- 
ably, thus the right upper corner of the table 2 is not physical. We can conclude 
that the formula ( p9| ) approximates the equilibrium concentration of monomers 
with relative errors less than 1%. 
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4 Solution of the parametric equations for Grinin's 
model of spherical micelles. 

There is another model of spherical micelles elaborated by Grinin ||T^. The 
hydrophobic tails of amphiphiles are assumed to be enough ordered in the core 
of micelles (on the contrary, for the drop model we assumed that hydrophobic 
tails are completely flexible as a liquid). The geometric restrains imply that the 
hydrophobic contribution to chemical potential of one amphiphile in aggregate 
is proportional to j/^/^. The hydrophilic effects can be described by the term 
proportional to v (see [0). The more detailed description can be found in the 
original article fl^ . 

For our purposes the final expression for the aggregation work is important, 

G. = ^(.^-l)-y(.-^/^-l) + g(.-l), (31) 

the two parameters^ are : a\ = a and a-i — b. The expression (|3l] ) is the 
decomposition of the aggregation work by the powers of v with p = 1/2 and 
with the coefficients 



2a b a? a? 2a b 

'° = y-2-4fe' = '' = Tb^ '" = -Y' ''=2- 

After calculating matrix of derivatives Cm,k and substituting it in (pO|), we obtain 
two equations. 



J o 3/2 



da 3 2bvs - avl'^ ' 



,2 



dvs „ 

db 2bv, - avV^ 



(32) 
(33) 



(in the nominator we have neglected a small term [a /2b — 2/3) by comparison 
with i/^^'^ and vl). The substitution 

Vs^{a/bf^\a,b) (34) 
reduces these equations to the quadratures 

2da _ dip{2ip - 1) 3db _ dip{2ip - 1) 
^ ~ tp{ip~2/3)' '2b ~ ip{ip - 2/3)' 



As for the drop model, the parameters a and b are not physical, but they are connected 
with real physical characteristics, 

where B is empiric constant (which equals to 1.4 under normal conditions, see Q), d is 
minimal distance between hydrophobic tail and water molecule, du^o characteristic size 
of water molecule, 1^ is hydrocarbon length, s is effective surface area of hydrophobic tail, 
<5 is characteristic scale of space distribution of charge formed by aggregate, ez is effective 
charge of hydrophilic head of amphiphile, sq is dielectric constant, e is dielectric permittivity 
of solution, ric is number of carbon atoms in hydrophobic tail. 
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which give after integration two equations of fourth degree connecting Lp with a 
and 6, 



^3|^-2/3| = 72(a)&', 

where 7i(&), 72(a) are arbitrary functions. Equahng the right-hand sides, one 
obtains : 71(6) = Cfo"^, 72(a) = Ca^'*, i.e. the equation is 

ip^\ip-2m^Ca-%^, (35) 

with an arbitrary constant C . Using small variations from the initial parameters, 
one can conclude that the factor ip^ changes a littlest Then we obtain an 
approximate solution 

95 = 2/3 + {CalhQ^)a-%^, (36) 
or, with the help of (jsj), we have 

v,^{2a/ih + ^G{N)a-^h'f, (37) 

where a constant 7^ can be defined for any N , and once defined (by the nu- 
merical solution of initial equation (||) for certain fixed parameters oq, 60)1 it 
becomes universal characteristic for the model under consideration. In the table 
3 we show the values of ^g{N) for some N . Now one can use the equation (||) 
to calculate ni. 



ni = exp 



^(l+7G(A^)a-^6^)^ 



(38) 



Thus we have found the complete dependence of the equilibrium concentration 
of monomers on the parameters for Grinin's model of spherical micelles. 

Note that the characteristic values of the parameters a and b of Grinin's 
model of spherical micelles are 

a«1.5, 6«0.1. (39) 

With the help of numerical simulations, we calculate the values of the constant 
7g for some N for these characteristic values of the parameters a and 6, see 
table 3. 



^® Indeed, using an explicit expression ( |3l[ ) for function Gu we obtain an equation for 
extrema, 

bv - a!/^/^ -I- /46 = In m , 
with larger solution (corresponding to the minimum Vs) 

Us = (a/26+ ^(lnnl)/6)^ 
whence it follows that ip > 1/2, therefore all changes occur in the vicinity of the point 2/3. 
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N 


2.0 


2.5 


3.0 


3.5 


4.0 


4.5 


5.0 


m 


1.6280 


1.6460 


1.6549 


1.6608 


1.6652 


1.6687 


1.6716 




-592.08 


-541.82 


-517.62 


-501.66 


-489.76 


-480.29 


-472.42 



Table 3. Grinins model: values of the constant jg CL^d of the equilibrium con- 
centration of monomers ni for some N calculated with the parameters a = 1.5, 
b = 0.1. 

The table 4 contain^^ the relative errors (in per cent) of the expression ( p8| ) . 



b \ a 


1.0 


1.1 


1.2 


1.3 


1.4 


1.5 


1.6 


1.7 


1.8 


1.9 


2.0 


0.05 


0.0 


0.9 




















0.06 


0.7 


0.0 


0.8 


















0.07 


1.2 


0.8 


0.1 


0.8 
















0.08 


0.9 


1.3 


0.9 


0.1 


0.8 














0.09 


1.3 


1.2 


1.4 


0.9 


0.0 


0.9 












0.10 




0.2 


1.4 


1.4 


0.8 





1.0 










0.11 






0.6 


1.6 


1.4 


0.7 


0.1 


1.2 








0.12 








1.1 


1.5 


1.3 


0.6 


0.3 


1.5 






0.13 








0.4 


1.5 


1.7 


1.3 


0.5 


0.5 


1.8 




0.14 










0.5 


1.7 


1.7 


1.1 


0.3 


0.7 


2.2 


0.15 












1.1 


1.8 


1.6 


1.0 


0.0 


1.0 



Table 4- Relative errors (in per cent) of the expression (|gqj calculated for a 
wide domain of the parameters a and b (N — 2). 

The absence of the values in the left low corner is caused by the condition 
ni > 1 (otherwise there are no stable aggregates). The right upper corner is 
empty by virtue of the condition ni < N (the opposite case has no physical 
sense). We see that the relative error of the approximate solution (^8|) changes 
in an irregular way, however, it does not exceed 2% when the parameters a and b 
take all physically reasonable values. This accuracy can be considered as a high 
one, thus we conclude that the formula (|3^) is a good approximation for the 
equilibrium concentration of monomers in Grinin's model. Nevertheless, there 
are problems which demand more accurate results. For this reason we would 
like to stress that the form of our expressions is explicit, therefore it is possible 
to analyze the results qualitatively and to compare different models. 

Note that the form of the solution ( |38| ) is rather different from drop's one 
(p9|). However, we can not compare them directly because the parameters a, 
b and 6i, 63 are completely incompatibles. Thus, one should express these 
parameters in terms of the physical quantities (such as the global characteristics 
of the aggregation work), substitute their into ( p9| ) and (p8|), and after that to 
compare the results. 



^"^ To fill this table, the equilibrium concentration of monomers ni for some values of 
parameters a and b was calculated by two different methods : by numerical solution of the 
equation and with the help of the approximate formula . 
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5 Conclusions 



Based on the nucleation theory, we have developed the method which allows to 
obtain the parametric equations under rather general assumptions : 

- closure of the system ; 

- dilute solution of surfactants ; 

- general model expression of G^, as a sum of the fractional powers of v ; 

- comparability of the numbers of amphiphiles in monomers and micelles. 
To deduce the parametric equations we have used these assumptions, but 

they are of the rather different significance. More important hypothesis is the 
closure of system. Some other ones are necessary to be mathematically rigorous. 
Finally, there are certain conditions which are used to specify the differential 
equations, to obtain concrete results (such as the solutions ( p9| ) or (^)), but 
they are not strictly necessary for the method's work. For example, the decom- 
position of function Gi, by the fractional powers of v is physically reasonable and 
it simplifies the method significantly. But this assumption can be diminished, 
i.e. the expression of Gi, can be generalized. In other words, by concretizing 
the form of some expressions (for example, (||), (|4|), etc.), we simplify mathe- 
matically our problem, physically leaving it in general form. 

We have obtained the differential equation for the important characteristic 
of micellization ~ for the aggregation number of micelles. Taking into account 
the remark in second part about the invariance of the parametric equations by 
addition of some function g{i/) to the aggregation work, we can conclude that 
Vs is an universal quantity of micellization. On the contrary, the equilibrium 
concentrations of monomers and micelles essentially depend on the concrete 
expression G^, consequently, they are not universal. 

The system of parametric equations was solved for two practically important 
model of spherical micelles : the drop and Grinin's models. We have obtained 
an explicit dependence of the equilibrium concentration of monomers on the 
parameters of micellization. The numerical simulations showed that these for- 
mulae have a high accuracy. 

At last, there appears a real possibility to verify the hypothesis that the 
physically measured quantities of micellization (e.g., the relaxation times) are 
independent of the concrete analytical expression of the aggregation work, but 
they are defined by its global characteristic (such as point of minimum, height 
of activation barrier, etc). If this hypothesis will be confirmed, it will justify 
the use of the simple model expressions of the aggregation work (such as ( p2| ) 
for the drop model, or ( |3l| ) for Grinin's model). Otherwise, our method allows 
to compare different model expressions and the solutions based on them, and 
to determine the substantiality of the model and its accuracy. 

We can conclude that elaborated method is effective enough both for applica- 
tions (for example, to determine the dependencies of equilibrium concentration 
of monomers on the parameters for the two models of spherical micelles) and 
for the general analysis. We hope that it will allow to solve some problems in 
the theory of micelles, and to comprehend profoundly the mechanism of micel- 
lization. 
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